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Abstract 

The closed system of equations for a determination of parameters distribution 
function for random interaction fields was calculated analytically. The estimations 
of critical concentrations for phase transitions in diluted crystal magnets with face 
centered cubic (fee), volume centered cubic (vec) and simple cubic (sc) lattices with 
interaction of Ruderman-Kittel-Kasua-Yosida between spins were made. 

Magnetic properties of disordered systems are investigated for a long time. In this 
direction great number of works were made, there are large number of books [1-5] and 
reviewed publications [6,7]. The most of approaches to theoretical investigation of mag- 
netism of disordered matter use a assumption about random distribution of exchange 
integrals in Hamiltonian of spin-spin interaction, at that a parameters of distribution 
function must have agreement with experimental data. We have formulated in [8,9] sev- 
eral different approach, which one help to calculate a distribution function of random 
fields of exchange interaction for amorphus magnets. A parameters of such function are 



in agree with each other and are calculated with using of law of exchange interaction (or 
magnetic moments of particles, clusters, grains etc.). In this article the proposed approach 
was extend for investigation of magnetic properties of crystal materials. 

1. DENSITY OF FIELD DISTRIBUTION 



In [8] the density of distribution of random field exchange interaction have been calculated 
for amorphus ferromagnetic material. In this case the probability of the hit of particle 
in element with volume dV can be determinate as and density of distribution of 
particles over magnetic moment S 
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For crystal lattices difference is in that the density of distribution of interacting particles, 



which ones have same magnetic moment So, over in the model of Ising is 
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Here a fc - - a relative number of particles, which ones have orientation in "positive" 
direction, f3 k — "negative", a k + k = 1, S(flk) — delta-function of Dirac, N — number of 
sizes of crystal lattice, N — number of interacting particles. The characteristic function 
in this case is 

A(p) = [ W{H) exp (ipH) dH = 
= 11 K 1 - P) +P («exp {dpp k {S , r M )j + exp {-ip(p k (S , r M )})] , (3) 
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where p = — , tp k (So,r k>0 ) — field, which one was created in origin of coordinates by 
particle in cite r fc - If to leave only three terms of expansion of exponent then 

In A{p) ^i{a- 0)pp ~ ^P [1 + (« - Pfp] p 2 ^ ~ 
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The distribution function for random field exchange interaction W(H) (it is same as in 



case of amorphus magnet) is "extended" (^-function with view 
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Remind that for amorphus magnets 



H — n I ipdV, 
v 



where n — volume concentration. 



B 2 « 2n / (fdV, 
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After thermodynamical averaging of characteristics of field "origin" in cite 
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|S| = So 
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and after configuration averaging for relative magnetization M per one cite, for < H > 
and < B 2 > easy to receive a system of selfconsistented equations 
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This system is analog of Sherrington-Kirkpatrick's system of equations, see for example 
[10,11]. Here B — parameter of order (similar to q), and in case of M — it describes a 
spin glass state. The advantage of system (8) is in that the basic parameters of distribution 
function (< H >, < B 2 >) are connected with each other and depend of law of interaction 



of particles. Here and future for < H > and < B 2 > we shall have left out a sing of 
averaging. 

The system (8) can be essentially simplified, if to make a substitute of distribution 
function 

by rectangular 

_ 0, -B > x, B < x 

W(x) = I (10) 

There are example of numerical solving of equation for M with exact and approximated 
functions [9]. Near of points of phase transition (small M and B) error in calculations 
is negligible. The available estimations of critical density p c are addition argument in 
favour of this substitution. p c corresponds to percolation threshold. It can be received if 
to consider of system (8) for case of direct exchange. 

So, from equation 
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it follows that a Curie point can be determined from relation 

§tanh(^H = l. (12) 
B \ kT J v ; 

It is obviously that condition — = 1 defines a maximal concentration p c , below of which 

B 

a ordering does not exist even at T — 0. For direct exchange <pk = fo, and summation 
must goes over nearest neighbors. 

= ^0 _ Peffo = l 
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where z — number of nearest neighbors. 

The critical concentrations p c , calculated by help of formula (13), have good agreement 
with knowing results in theory of percolation [12-15]. 



2. RKKY INTERACTION 



For RKKY interaction the "strength of exchange field" is 

<P = AF(x), (14) 

where x = 2kpR, kp — impulse of free electron on Fermi surface, 

. xcosx — sinx 

F ( X ) = 14 ' ( 15 ) 

and A has dimension of magnetic field and defines of intensity of exchange interaction. 

Since for "standard" metal kp = 3n 2 n s , where n s — concentration of free electrons, 
Fermi-impulse is comparable with parameter of lattice (kp ~ a). The results of summation 
in formulas for Hq and B 2 can have significant dependence from mutual disposition of 
interacting atoms, i.e. from type of crystal lattice. It is known that distance between 
atoms for different lattices can be calculated as 

1) SC 

R ni ,n 2 ,n 3 = ay/nl + n% + nl, (16) 



2) FCC 



2) VCC 



R ni ,n 2 ,n 3 = | + n 2 ) 2 + (ni + n 3 ) 2 + (n 2 + n 3 ) 2 , (17) 



Rn 1: n2,n 3 = ( n i + n 2 ~ n 3 ) 2 + Jn-y + n 3 - n 2 ) 2 + (n 2 + n 3 - n x ) 2 , (18) 

where ni, n 2 and n 3 — whole numbers. The number of atoms, which ones have same 
distances, equals to number of integer roots of equations (16-18). 

The parameters of distribution function Hq and B 2 for crystal magnetic alloys, which 
ones have given type of lattice, can be calculated by summation: 

H =pA F(2k F R nun2 , n3 ), (19) 

n\,ni,n 3 
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B 2 = 2pA 2 F 2 (2k F R nun2 , n3 ). (20) 

ni,n 2 ,n 3 

At the calculation of H and B 2 it is necessary to have in view, that concentration of free 

4 2 1 

electrons n s = for FCC lattice, for VCC n s = — and for SC lattice n s = — . 

a 6 cr or 

The numerical estimations of parameters of distribution function W(H), which ones 
have been made with taking in account a first 20 roots of equations (16-18), when sums 
(19-20) go to suturation, for FCC, VCC and SC lattice, correspondingly, are: 

TT 

FCC: H = 0.037pA, B = 0.012^^4, = 3.14^ ; 

B 

TT 

VCC: H = 0.032pA, B = 0.011 JpA, = 2.83^ ; 

B 

TT 

SC: Hq = 0.011pA, B = 0.016^.4, = 0.69 ^/p . 

B 

From this data it follows that for FCC lattice in case of RKKY exchange in standard 
metal the ferromagnetic state is possible at p > p c = 0.32, and the spin glass ordering 
takes place at p < p c . For VCC the ferromagnetism is possible only at p > p c = 0.35. And 
for SC it is possible only the spin glass state at low temperature and the paramagnetism 
at high one. 

Thus, the method of random field exchange interaction allows to make a estimations 
of critical concentrations of phase transitions in diluted magnets in dependence of crystal 
structure and of type of exchange interaction. 
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